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Abstract The paper deals with the explicit calculus and the properties of 

the fundamental solution _ftr of a parabolic operator related to a semilinear 

CN| , equation that models reaction diffusion systems with excitable kinetics. The 

initial value problem in all of the space is analyzed together with continuous 

dependence and a priori estimates of the solution. These estimates show that 

Q,' the asymptotic behavior is determined by the reaction mechanism. Moreover 

I I it's possible a rigorous singular perturbation analysis for discussing travelling 

1^ . waves with their characteristic times. 
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■"nI" ' 1 Statement of the problem and results 

•^1 Let f2T = {{x,t) : X e ^ < t < T} and, for functions u{x,t) € C^{f2T ), 

CO . let Vq be the non linear initial value problem defined in all of the space: 

o 

Lu ^ ut — euxx + au + b L e~"^*^'^^ u{x, t) dr = F{x, t, u) 



u{x,0) = g{x) a; e 5R, 



(1) 
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where a, b, e, /3 are positive constants and F, g are known functions of their 
arguments. 

If K{x, t) is a fundamental solution of the parabolic operator L and 
F(x, t, u) verifies appropriate assumptions, then the differential problem 
(H)) is equivalent to the integral equation 

u{x,t) = F[u{x,t)], (2) 

where ^ [ w ] is the mapping 

T[v{x,t)] = / K{x-^,t) g{Odi + (3) 

+ [ dT f Kix-^,t-T)F[C,T,vii,T)]d^. 

Jo Ju 
Now, let ||w||t = sup |ti(a;,i) I, and let Bt denote the Banach space 

Bt = {vix,t) : v £C{Qt), \\v\\t < oo }. (4) 

The aim of the paper is the explicit calculus of K and the analysis of its 
basic properties (§2) useful to prove that J^ is a contraction of Bt in Bt 
and hence admits a unique fixed point u — u{x,t) € Bt- Thus existence 
and uniqueness results for the problem Vq are deduced (§3), together with 
continuous dependence and a priori estimates of the solution (§4). 

The following estimate (§2) is worthy of remark 

— -^ —at —Bt 

\K{xM < 7rr=, [e~"* + bt —. ] (5) 

Z V 7r£t p ~ a 

and it shows that K has the same basic properties of the fundamental so- 
lution of the heat operator and moreover decays exponentially to zero as t 
increases. As consequence (§ 4) , the solution u of Vq is such that 

Mx,t)\\T < \\g\\ (1 + TrVbt) e-"* + /3o ||F||, (6) 

where a; — min(a,/3), the norms are defined in §4 and the constant /3o 
depends on a , b, f3 (see ([50]) ). So, when t is large, the asymptotic behavior 
of u is determined by the properties of the source F. 

The equation ^ models several engineering applications as motions of 
viscoelastic fluids or solids [51I^ [T^[^ : heat conduction at low temperature 
[TSII^ . sound propagation in viscous gases [TB] and perturbed sine Gordon 
equation in the theory of the superconductivity [2l[7l[28l or in the propagation 
of localized magnetohydrodinamic models in plasma physics [.30] . 
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As an example, in section 5 the results are applied to the FitzHugh - 
Nagumo equations (FHN) which model many important biological phenom- 
ena pi fTTl[Hl[T51[^ ■ According to well known results concerning systems of 
non linear reaction- diffusion equations P^I5T| . estimates like ([5]) confirm 
that the large time behaviour of the solution is determined by the reaction 
mechanism. Moreover by means of the explicit fundamental solution K and 
its properties, it's possible a rigorous singular perturbation analysis to ap- 
proximate travelling pulses in excitable media, together with their transition 
times. 



2 E\indamental solution and properties 

By denoting with 



u{x,s) = / e u{x,t)dt , F{x,s) = e F [x,t,u{x,t)]dt , 

Jo Ja 

the Laplace transform with respect to t, from ^ we get 



u{x,s)= / Kix~C,s)[9ia + Fi^,s)]dC, (7) 



where 



K(x,s) — r= — with a = s + a + . (8) 

2 y/e a s + p 

Therefore, if K{x,t) represents the inverse Ct transforms of K{xts), 
from ([7]) formally it follows that 



u{x,t)^ / K{x-^,t)g{i)d^ + (9) 



+ ( dT f Kix~C,t~T)F[^,T,ui^,T)]dC 
Jo J^i 

If r ~ \x\ /y^ and J„(z) denotes the Bessel function of first kind and 
order n, let us consider the function 



Kir,t) ^ -^= e-°'- (10) 

2y'7:et 



e-Pi*-v)j,{2^byit-y) ) dy. 



2 V 7re 7o \/i-y 
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Theorem 2.1- In the half-plane 5ftes > max{~a, — /?) the Laplace in- 
tegral Ct K{r, t) converges absolutely for all r > 0, and it results: 



CtK= e-'' K{r,t)dt ^ —^. (11) 

Jo ^yjSCT 

Proof - For all real z, one has | J„ (z)| < 1 and the Fubini -Tonelli 
theorem implies that 



LtK = 



2^e{s + a) 



^ ^°°e-(^+")^-5|dy r e~{s+0)t J,(2y/b^t) '^^ 



2%/7re 7o Jo ^ft 

Since 

/ e-P* ,/^ Ji(2V^) dt = 1 - e-'^/P (SRep>0), 

it follows that 



K (r, s) = - — -= / e *y ^ ^ ^ s+ft i » 



2 \/7r£ 7o V^ 2 7e cr 



Theorem 2.2 - T/ie function K has the same basic properties of the 
fundamental solution of the heat equation, that is : 

A- K{x,t) e C°° for t>0, X e 3?. 

B - For fixed t > 0, K and its derivatives are vanishing esponentially 
fast as \x\ tends to infinity. 

C - For any fixed d > 0, uniformly for all \x\ > S, it results: 

lim K{x,t) = 0, (12) 

D. - For t > 0, it is LK ^ 0. 

Proof - By (fTU)) the properties A-B-C are obvious. To verify D we put 
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by 



^{y,t) = Jj^ M2^by{t-y) )e-^(*^«) (14) 

Kir,t) = V(r,i) - / i^{r,y) ip{y,t)dy. (15) 

As 

{dt+ a- drr) i^{r,t) = (16) 

from (ITSl) it results: 



Kt+aK ~ Krr = ~f{t, t) V' (r, t) - j [ip{^t + aip) - (pi/jrr]dy = 

Jo 

= -(p{t,t)^p{r,t)- [^{(fit+aip)- ip{tpy + aip)]dy. (17) 

Jo 

But 

t pt 

ipipydy ^ip {t, t) i;{r,t) - ip ipy dy 



and so (fTT)) gives 

Kt + aK-Krr = - ^ (ipt + iPy) dy. (18) 

Jo 

As 



{dt+dy)^{y,t) = b M2^by{t-y)) e-'^'*"^) (19) 

from ([T5)) it follows 

JTt + aX - X„. = - 6 iiTi (20) 

where Ki is given by 

K,{r,t) = -^ r e-5|-s.-/J(*-y) J,(2 ^by {t ^ y)) ^. (21) 

2 0^ Jo Vy 

On the other side, the convolution e~'^* * iiT is 



'0 

and moreover it results: 



K ^e-^UiJj - / 'ilj{r,y)dy / e' '^'^^'^'^ Lp{y,T) dr (22) 
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by 



^(*-^)^(y,T)dT = e-''(*-^) / J^ J,{2^by{T^y))dr 



T-y 



(23) 



-/3(t-j/) 



1 - Joi2^byit-y)) 



As consequence, from ((22|) - ((23l) we get 



e-l'UK^ e-^ (*-a) ^ (r, y) Jo ( 2 v/6 y (i - y) ) dy = K^ (24) 

Jo 

and so (|20p - ([M)) imply property D because i^T^r = £ Kxx ■ B 

Moreover, if /„(z) denotes the modified Bessel function of the real vari- 
able z and 



E{t) 



>0. 



(25) 



a - /3 
then the following theorem holds: 

Theorem 2.3 - The fundamental solution K (r, t) satisfies the estimates 



\K\ < 



2V7ret 



btE{t)] 



(26) 



f \K{x~(,t)\d^ < e-'^' + VbTTte-^* Io(^^t\ (27) 



Proof - As \Ji{2 y/by{t - y) \ < y/by{t~y), from ^ we get 
Further, because 



2Vnet 



d^ - 1, 



from ([TU)) and |Ji(^:)| < 1 it results: 

/ \K{x-^,t)\d^ < e--' + [ 
J?R Jo 



-ay-Pit-y) ^^y 

\ t-y 



+ {^/b/2)Trte^^ 



Io{'-^t]-lJ^-lt 
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As Ii{-z) = -h{z), Io{-z) = Io{z) and Ii{\z\) < Io{\z\), the estimate (p7|) 
follows. I 



Now, if one puts w = min (a, /3), by means of (PT)) we get 

/ |if(a;-^,i)|d^ < e-"* + \/67rte-'^* (28) 

because /o(|-z|) < ea;p(|z|). As consequence, for all t < T, from (P5| it follows 
that 

[ dr [ \K{x-£.,t-T)\d£. < {l+VbTr/uj)T (29) 

while, for i — > oo, the formula (|27p implies 

f dr f \K{x-C,t)\d^ < - + TrVb f re"^^ ^o(^^^t) dr 
Jo Jk o, Jo 2 



i+.y6 « + ^ 



a 2(a/3)3/2 



(30) 

= f3o. 



The integrals related to the kernel Ki = e ^* * K defined by (pij) can 
be estimated in the same way. As | Jo| < 1 and 

r E{T)dT = — =p^, (31) 

Jo "'P 

the following results can be stated: 

Theorem 2.4 - For a/Z t £ (0,oo), the kernels K, Ki defined by (fTOl) . 
(|2ip are smc/i t/iat 

[ \K\di< {1 + VbTTt)e-''\ [\Ki\di<E{t) (32) 

Jk J5R 



ft !• rt 



f dr f \K\dC< /3o, f dT f \Ki\ d^ < A (33) 

Jo JSR Jo Js? 

where uj = min{a, /3), E{t) is given by (j25p and the constants /3o, /3i are 

defined in §U^, §1^. I 
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At last we observe that 



Ki{x - 1 1) d^ = x(i), / K{x - e, t) d^ = (dt + (3) x{t) (34) 



with the function x(i) defined by 



i+± f , s , 1 



X{t) = e- — ' senigt)/g, g^ - ^46 - (a - /3)2. (35) 

The formulae (I34p . psp can be verified by means of their C transforms be- 
cause the related Laplace integrals are endowed with absolute convergence. 



3 Existence and uniqueness results 

As for the data F and g of the problem Vq , we shall admit: 

Assumption 3.1 The function g{x) is continuously differentiable and 
bounded together with g'{x). The function F{x,t,u) is defined and contin- 
uous on the set 

D = {{x, t, u) : [x, t) e Qt , ~oo < u < oo} (36) 

and more is uniformly Lipschitz continuous in {x, t, u) for each compact 
subset of Qt ■ Moreover F is bounded for bounded u and there exists a 
constant Cp such that the estimate 

\F{x,t,ui) - F{x,t,u2)\ < Cf |ui-U2| (37) 

holds for all [ui, U2). M 

When the problem Vq admits a solution u, then u must satisfy the 
integral equation 



u^ f Kix~t t)giOd^ + f dr f Kix-C,t- T)F[t r, w(f , r M (38) 
J?R Jo JU 

because of the properties of instated by theorem 2.2 and the assumption 
3.1 on the data . On the other hand, let u{x,t) be a solution of (1551) which 
is continuous and bounded so that also F(x,t,u{x,t)) is continuous and 
bounded by assumption 3.1 . Then one can verify that u satisfies (jlj owing 
to the properties oi K , F and g. So, it' s possible to conclude that 
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Theorem 3.1 - The differential problem ^ admits a unique solution 
only if the integral equation (|38[) has a unique solution which is continuous 
and hounded in Qt ■ B 

Consider now the integral equation (1551) and, for small time < t < 9 
(with e < T), let 

||u||e = sup \v{x,t) I (39) 

-3c < :c < oo 
< t < 9 

and 

Bg = {v{x,t) : V e C [(~cx),cx)) x [0,6*]] and \\v\\g < oo } (40) 

Further, let denote by /i * /2 the convolution 



/l(x,i) */2(x,t) - / /i(^,t) /2(x-^,t)d^ (41) 



and by Tv the mapping 



Fv{x,t) = K{x,t) *g{x)+ K{x,t-T) ^F{x,T,v{x,T))dT (42) 

The foregoing remarks imply that Bq is a Banach space and (HU maps 
Be in Be and represents a continuous map for < i < 6* . In fact, from ([57|) 
and (1^ it results 

ll^i'i - -Ft>2)||e < Cf 6'(1 + 7r\/6/a;) Hvi-'yalle. (43) 

Hence, when one selects 6 such that Cp ^ ( 1 + 7r-\/6/w ) < 1 then J^ is a 
contraction of Be into Be and so has a unique fixed point u{x,t) e Be ■ 

In order to show that this result holds also for < i < T, it suffices to 
proceed by induction and so the integral equation (1551) has a unique solution 
for < t < (n + 1) 9, whatever the positive integer n may be. 

Theorem 3.2 - When the data {F, g) satisfy the assumption 3.1, then 
the initial value problem Vq defined in ([T]) admits a unique regular solution 
u{x,t) in [2t ■ B 



4 A priori estimates and continuous dependence 

According to the assumption 3.1, in the class of bounded solutions, let 



II5II =sup|5(a;)|, ||u||j, =sup|u(a;,i)|, ||i^ || = sup |F(x,i,M)| 
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The solution u(x,t) of the problem Vo depends continuosly upon the data, 
owing to the following theorem. 

Theorem 4.1 - Let ui, U2 be solutions of the problem Vq related to 
the data {gi, Fi) and {g2, F2) which satisfy the assumption 3.1. Then there 
exists a positive constant C such that 

\\ui -U2 llr < C sup |gi -52I + C' sup \Fi{x,t,u) - F2{x,t,u)\ , 
3? D 

where C depends on Cp, T and the parameters a, b, /3. 

Proof - The integral equations ([55)) for ui and U2 are differenced and 
the estimates of theorem 2.4 allow to apply the Gronwall lemma. I 

Remark 4.1 - When the function F is a known linear source f{x,t) 
then ((551) gives the explicit solution of the linear problem Vq 

u^K*g+ / K{x,t-T) =¥ f{x,T)dT. (44) 

"'0 

By means of this formula and the basic properties of the kernel K stated in 
theorems 2.3 and 2.4, it's possible to estimate u and its derivatives together 
with their asymptotic behavior. I 

In the non linear case, the integral equation psp implies a priori estimates. 
For example: 

Theorem 4.2 - When the data {g, F) of the non linear problem Vo 
verify the assumption 3.1, then it results: 

\Hx,t)\\^ < Pq \\F\\ + \\g\\ [e-°* +\/6^te~"*] (45) 

where the constant j3q is defined by pOp and depends on a, b^ /3 . 

Proof -It suffices to apply to (p8l) the estimates (|28)) and (l33l) . I 



5 Excitable models and travelling pulses 

As it is well known, a biologcal system is excitable if a stimulus of sufficient 
size can initiate a travelling pulse which will propagate through the medium. 
One of the most revelant examples of systems with excitable behavior is given 
by neural communications by nerve cells via electrical signalling |111I15| . 

Let u{x,t) a transmembrane potential and let v{x,t) a variable associ- 
ated with the contributions to the membrane current from sodium, potassium 
and other ions. Then a simple example of excitable models is the Fitzhugh 
Nagumo system (FHN) [2n[22] : 
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' du d^u 

(46) 

— - ^ bu - Pv 
y ot 

where £ > is the diffusion coifEcient related to the axial current in the 
axon, while h and /3 are positive constants that characterize the model's 
kinetic. Further 

f{u)^u{a-u){u~l) (47) 

In absence of diffusion (e — ), the system (l46l) - (l47)) is related to a two 
variable phase plane system with the phase portrait that varies according to 
different values of the parameters a, b, P . In fact, when (1 — a)^ < 46/^, 
there is only the stable steady state O with u = 0, w = 0. If (1 — a)^ > 
46//3, there are two other steady state A and B with 



UA = ■;^ [a + 1 - \/(l - a)2 - 46//3] < a, va ^ -ua 



and 



UB ^ -[a + 1+ A/(l-a)2 - 46//3] < 1, ^s = -^ us- 



The state B is as locally stable as O while the state A between them is 
unstable. So, a sufRciently strong perturbation to a point ( u, 0) with u > 
UA may induce a large phase traiectory excursion, while a small perturbation 
around O {u < ua) rapidly decays towards O . All this is typical of the 
threshold behavior. 

Moreover, ii z = x — ct, travelling wave solutions of (j46p are functions 
m(z),z)(z) such that 



eil + cii + f{u) -V =0, cv +bu - l3v =0. (48) 

A solitary pulse is meant to be a solution of (j48p with appropriate bound- 
ary conditions for z — ?> ±oo. A classical meaningful example [lUll^l^ is 
related to the particular case oi b — and 



lim u{z) = 1, lim u{z) — lini v{z) — 0. (49) 

In this case the exact analytical solution of (|48)) . (|49)) is 
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where: 7 = l/-\/2F, c = Y^e/2 (1 — 2a). Several other solutions of the 
equation (051)1 which behave like ([SU]) are well known in literature [HOIHl 
[T51[^[^ and are all bounded solutions that verify the assumption 3.1. As 
consequence, the results stated in section 3 and 4 are meaningful for the 
analysis of the evolution of (FHN) model in all of the space. 

Let 

u(a;, 0) = wq : v{x,0) = vq (x gSR) (51) 

the initial conditions related to the system (|iS)) and let 

f{u) ~ — au + (p{u) with ip = u {a + 1 — u). (52) 

As (|46l) 9 implies 

V = voe-^' + b f e^^^'-^'>u{x,T)dT, (53) 



when one puts F{x, t, u) ~ ^{u) — vo{x) e ^* , then the initial value prob- 
lem (|46)) - ((5lT) can be given the form ([T]) to obtain u in terms of the data. 

For this we apply the formula (pS)) with g — uq and put 

K(^F= [ dr [ K{x-^,t-T) F[^,T,u{^,T)]d^. 
Jo J?Si 

Because of ([Ml) it is iiT (g) (ug e ~ '^ * ) = vq* Ki, and so 



K®F = K®ip- Ki* vo{x). 
Therefore by (l38l) we get: 



u{x, t) ^uo*K — vo*Ki+ip(>^K (54) 

and this formula, together with (1531) . allows to obtain also v{x,t) in terms 
of the data. If we observe that 



Title Suppressed Due to Excessive Length 13 



t 
e-/3(*-^)Xi(x,T)dT = K2{x,t) 



I 



t f.-f^-'^y-Pii-y) jt-y 



by means of ((53)l - ((55|l we get 



y^ U2^by{t~y) )dy (55) 



v{x,t) ^ v^e ^ * + b[uo*Ki ~ vo*K2 + ^ ® Ki]. (56) 

When the non linear part •^{u) of f{u) is neglegible, or can be approximated 
by a constant (^o, then ((54|) - ((56)) represent the explicit solution (u,w) of the 
linear case. Otherwise (j54p is an integral equation for the unknown u that 
is identical to p8l) and then ([55)1 can be applied for v. 

Because of theorems 3.2 and 4.1, the following result can be stated. 

Theorem 5.1 - When the data {uq,vq) satisfy the hypotheses 3.1, then 
the initial value problem related to (FHN) system (I46p has a unique regular 
solution in the space of bounded solution. Moreover this solution depends 
continuously upon the data. I 

In order to apply theorem 4.2, an estimate like ([5^ must be set also for 
the kernel K2. From (1551) it results: 



[ \K2ix - ^,t)\ d^ < [ e-'^y-f'^'-yHt~y)dy <tE{t) (57) 

SO that, if 

||uo|| =sup|uo(2;)| , ||wo|| =sup|?;o(2:)| , \\lp\\ = s\Yp\ip{u)\ , 



by means of ([57)) and theorem 2.4 the formulae (|M)) - (1551) imply the following 
conclusion. 

Theorem 5.2 - For initial data {uo,vq) compatible with assumption 3.1, 
the solution of the problem (|46p - (j5ip satisfies the estimates: 

\u\ <\\u4 (l + 7rV6t)e-"* + ||«o|| E{t) + (i^ M\ 

(58) 
\v\ < Ill-oil e-^* +h{\\uo\\ +t\\v4)E{t) + b Pi M 

with Po, Pi, E{t) defined &?/ ([301), dSH), ([251). I 
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The analysis and the stabihty of solutions of nonlinear binary reaction - 
diffusion systems of PDE's, as well as the existence of global compact attrac- 
tors, have been discussed in a great number of recent and interesting papers 
(see e. g. [SJirz, 26,27,31]). As it is well known, the (FHN) system admits 
arbitrary large invariant rectangles E containing (0, 0) so that the solution 
{u,v), for all times t > 0, lies in the interior of S when the initial data 
{uotVo) belong to S . Besides, when both /? and a are strictly positive, 
then the estimates ([55]) confirm that the large time behaviour of the solution 
is determined only by the reaction mechanism because the explicit terms 
depending on the initial data are exponentially vanishing. 

Finally we observe that when 5 = 0, /3 = there is no bounded set that 
attracts all solutions and the asymptotic analysis needs the use of character- 
istic parameters such as fast and slow times, together with rigorous estimates 
uniformly valid for all t. To this extent, the basic properties of the kernel K 
might be useful for a rigorous singular perturbation analysis pTHfTHI^ . 
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